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We investigate the Mott effect for heavy quarkonia due to Debye screening of the heavy 
quark potential in a plasma of massless quarks and antiquarks. The influence of residual 
color correlation is investigated by coupling the light quark sector to a temporal gauge field 
driven by the Polyakov loop potential. This leads to an increase of the Mott dissociation 
temperatures for quarkonia states which stabilizes in particular the excited states, but has 
marginal effect on the ground states. 



1 Introduction 

Since the suggestion of J/ip suppression as a signal of quark-gluon plasma (QGP) formation by 
Matsui and Satz [T] in 1986 the problem of quarkonium dissociation in hot and dense strongly 
interacting matter has played a key role for QGP diagnostics in relativistic heavy-ion collision 
experiments. The original idea was that in a QGP the string tension of the confining potential 
vanishes and the residual one-gluon exchange interaction undergoes a Debye screening by the 
color charges of the plasma. When the temperature dependent Debye radius ro{T) (the inverse 
of the Debye mass m]j(T)) becomes shorter than the Bohr radius of the charmonium ground 
state (J/i/j) then the Mott effect [5] (bound state dissociation) occurs and the corresponding 
temperature is T^^^. This simple idea grew up to a multifacetted research direction when not 
only in the first light ion - nucleus collisions at the CERN NA38 experiment, but also in proton 
- nucleus collisions at Fermilab J/V-J suppression has been found so that there is not only a QGP 
but also a cold nuclear matter effect on charmonium production, see [5] for a recent review. 

If one wants to explore the question of screening in a plasma more in detail then a variety 
of appraoches is available in the literature, from the original Debye-Hiickel approach |3j where 
one can study any vacuum potential (for example the Cornell potential) and see its medium 
modification, over the thermodynamic Green functions approach to the ab-initio studies of 
heavy-quark potentials in lattice QCD. With the obtained medium-dependent potentials one 
can then study the bound state problem by solving the nonrelativistic Schrodinger equation or, 
more systematically, the thermodynamic T - matrix for quarkonia [5] . 

On the other hand one may calculate proper correlators directly from lattice QCD and 
extract from them spectral functions [6] . There is an intriguing disagreement between the Mott 
temperatures deduced from these spectral functions and those of the potential models: the 
latter are much smaller than the former! From the lattice data for quarkonium correlators one 
has extracted r^°" « l.gr^ while in potential model calculations T)J^" 1.2Tc. This problem 
has lead to the discussion of the proper thermodynamical function to be used as a potential 
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in the Schrodinger equation. Should it be the free energy or the internal energy? We will not 
follow this question in the present work, but refer to [3l [7] and references therein. 

Here we examine a simple model of screening as derived from one-loop calculations in thermal 
quantum field theory, and make a small extension of this result by putting the internal fermion 
lines in a constant temporal gauge field which mimics confining gluon dynamics (Polyakov- 
loop potential). In our approach the medium is made of plasma of massless quarks, described 
by the chirally symmetric phase of the Nambu-Jona-Lasinio (NJL) model of QCD Con- 
finement is implemented in the most simple way by coupling the system to the Poyakov loop 
variable - resulting in the so called Polakov loop NJL- model (or PNJL model). Recently it has 
been successfully used to reproduce lattice data [9] or to describe light meson physics at finite 
temperatures and densities [10| . 



2 Debye-screening in a PNJL quark plasma 

Given the self energy (polarization function) of a boson that mediates the interaction, the 
screened potential is given by a resummation of one-particle irreducible diagrams ("bubble" 
resummation = RPA) [11] 

VsM^V{q)/[l + F{0;q)/q'], (1) 

where we take V{q) = — |.g^/<z^, = |qp as the unscreened vacuum potential. The longitudinal 
gluon polarization function F(0;q) = — noo(0;q) in the finite T case can be calculated within 
perturbative thermal field theory where it takes the form 

^oo(^c^^; q) =Tg^ / 7^ Tr[7°5$(iL^„; p)7°5$(zu;„ ~iu;i;p- q)] , (2) 

where uJi = 2-kIT are the bosonic and Wn = (2n+ 1)t:T are the fermionic Matsubara frequencies 
of the imaginary-time formalism. The symbol Tr stands for traces in color, flavor and Dirac 
spaces. 5$ is the propagator of a massless fermion coupled to the homogeneous static gluon 
background field (^3 . Its inverse is given by [9l [TOj 

5^^(p;w„) = 7 • p + 7oici^„ - A3i,53 , (3) 
where Lp^ is related to the Polyakov loop variable defined by [S] 

<i>(T) = iTr,(e''5^-^'^-^) = ^(1 + 2cos(/3^3)) • 

The physics of $(T) is governed by the temperature-dependent Polyakov loop potential U{^), 
which is fitted to describe the lattice data for the pressure of the pure glue system |9]. After 
performing the color-, flavor- and Dirac traces and making the fermionic Matsubara summation, 
we obtain in the static, long wavelength limit 

noo(0;q) = — / dpp — = ^ — / dppU{p) = /($) = -mj^iT) . 

(4) 

where mD{T) is the Dcbye mass, the number of degrees of freedom is A'dof — Nc Nf = 6 and 
f<s>{p) is quark distribution function [lU]. The screened potential is thus 

VUq) = -47ra/[q'+ml{T)]. (5) 
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In comparison to the free fermion case [HI US] the coupling to the Polyakov loop variable 4>(T) 
gives rise to a modification of the Debye mass, given by the integral 







In the limit of deconfinement ($ = 1), the case of a massless quark gas is obtained (/(I) = 1) 
while for confinement ($ = 0) one finds that /(O) = 1/9. For the temperature dependence of 
$(r) we employ in the following chapter the results of a nonlocal PNJL model [13]. 



3 Variational ansatz and estimation of Mott temperatures 

Here we will use the derived potential in the quantum mechanical way to estimate the dissoci- 
ation temperature. With the trial wave function for the IS state 

V'is('';7) = y^exp(-7r) (7) 
and the non-relativistic two-body Hamiltonian 

rriQ r 

where the potential term is the Fourier transform of screened Coulomb potential ([5]), we obtain 
the energy functional for the Ritz variational principle 

E,,i,,T) = (^,s(7) I H I ^,s(7)) = ^ - ^^^'r)l2,y ' 

Simultaneously satisfying the conditions for the ground state energy dEis{j,T)/d'-f = and 
for a vanishing binding energy (Mott effect), Eisij,T^°*^) = 0, provides us with an ana- 
lytic expression for the critical Debye mass to_d(T'-J^°") = 27 . Using once again the condition 
^^is(7,T"°") = 0, results in the Mott condition for the Debye potential [2] 

r^(T^^°")=ao, (10) 

where ao = 2/{amQ) = l/^eaniQ is the Bohr radius and £0 — a'^mq/A the binding energy of 
ground state in the vacuum (mu =0). Due to the temperature dependence of the Debye mass, 
we obtain the Mott dissociation temperature in the massless quark gas (for <I> = 1) 



= ^3eo™Q/^dof/5 = ]l y^eomy (27r7Vd^ . (11) 

In Tab. |T] we give the parameter values according to set (i) of Ref. [H]. The Poyakov loop 
variable contribution considered in previous section affects the Mott temperature in Eq. (jlOp 
and gives the following formula 



nMott,* _ rriMott 



r"°"//F(¥) . (12) 
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Table 1: Mott temperatures T'^"" (yMott,*-) according to 
Eq. pT|) (Eq. p2)) ) for a massless ideal quark gas (PNJL 
model). The critical temperature is Tc = 202 MeV [l3| . 
The parameters are fixed to reproduce quarkonium states 
in vacuum as Coulombic bound states [14j . In the char- 
monium (bottomonium) system the heavy quark mass mg 
is rric = 1.94 GeV (mf, = 5.1 GeV) and the ground state 
binding energy eo is 0.78 GeV (0.75 GeV). 
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Figure 1: Temperature dependence of the Debye screening radius roiT) for charmonium (upper 
panel) and bottomonium (middle panel) with (solid lines) and without (dashed lines) coupling 
to the Polyakov loop $(T) shown in the lower panel (from [E]). Sequential dissociation of 
quarkonia states occurs at the Mott temperatures T^°^^ for which their Bohr radius equals 
roiT). 

This means that in the case of confining color correlations (0 < $ < 1) the Debye screening 
radius is larger than in a free quark gas at the same temperature, so that bound states get 
stabilized against thermal dissociation by color screening. 

The influence of the Poyakov loop variable on the dissociation temperature is summarized 
in Figure [1] which shows the temperature dependence of the Debye radius ro(T) = l/muiT) 
compared to the Bohr radii of the low-lying states of the charmonium and bottomonium family, 
respectively. Due to its larger mass and smaller Bohr radius, the T dissociates at higher 
temperatures than the J ftp, where the free quark gas case is almost reached and the Mott 
temperatures which almost coincide for the cases with and without coupling to the Polyakov 
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loop field. For the lighter J/ij} there is a noticeable stabilization due to the coupling to the 
Polyakov loop potential which results in an increase of the Mott dissociation temperature from 
1.25 Tc to 1.37 Tc; still more similar to results of nonrelativistic potential models rather than 
the still higher dissociation temperatures conjectured from the spectral functions deduced from 
lattice data for heavy quarkonium correlators by the maximum entropy method. For estimating 
the Mott temperatures of the excited quarkonia states we have employed the scaling of bound 
state radii: r^^ = 1.5 rj/^, r^t — 2.0 rj/^ as obtained in the Cornell-type potential model [J- 

4 Conclusions 

We have applied the methods of thermal field theory to estimate the effects of Debye screening 
on heavy quarkonia bound state formation. In order to account for residual effects of confining 
color correlations in the deconfined phase, we have used the PNJL model in the evaluation of the 
one-loop polarization function. As expected, a stabilization of bound states in the vicinity of the 
critical temperature for T > T^'is obtained. We applied Ritz' variational principle to derive the 
Mott criterion for bound states of the statically screened Debye potential and obtained Mott 
temperatures in good agreement with previous results from nonrelativistic potential models 
exploiting lattice QCD singlet free energies as potentials in the Schrodinger equation for heavy 
quarkonia. This agreement with previous results (see, e.g., Ref. [7]) includes also higher 
quarkonia resonances for which the stabilization effects is more pronounced. 
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